February 2, 2008 



The Effect of Cosmic Ray Diffusion on the Parker Instability 

Dongsu Ryu 1 , Jongsoo Kim 2 , Seung Soo Hong 3 and T. W. Jones 4 

ABSTRACT 

The Parker instability, which has been considered as a process governing the 
structure of the interstellar medium, is induced by the buoyancy of magnetic 
field and cosmic rays. In previous studies, while the magnetic field has been fully 
incorporated in the context of isothermal magnetohydrodynamics, cosmic rays 
have been normally treated with the simplifying assumption of infinite diffusion 
along magnetic field lines but no diffusion across them. The cosmic ray diffusion 
is, however, finite. In this work, we take into account fully the diffusion process 
of cosmic rays in a linear stability analysis of the Parker instability. Cosmic 
rays are described with the diffusion-convection equation. With realistic values 
of cosmic ray diffusion coefficients expected in the interstellar medium, we show 
that the result of previous studies with the simplifying assumption on cosmic ray 
diffusion applies well. Finiteness of parallel diffusion decreases the growth rate 
of the Parker instability, while the relatively smaller perpendicular diffusion has 
no significant effect. We discuss the implication of our result on the role of the 
Parker instability in the interstellar medium. 

Subject headings: cosmic rays — instabilities — ISM: magnetic fields — MHD 

1. Introduction 

In stability analyses of the interstellar medium (ISM), Parker (1966, 1967) put forward 
a simple model of the ISM which is composed of a single phase gas, magnetic field and cos- 
mic rays (CRs) under external, uniform, vertical gravity. It was assumed that gas pressure 
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originates from the ram motion of cloudlets, rather than the thermal motion of atoms or 
molecules, so the velocity dispersion of cloudlets was taken as the sound speed. In equilib- 
rium, the magnetic field has only a regular component and the ratios of pressures of magnetic 
field and CRs to gas pressure are constant. In addition, the CR dynamics was simplified 
by setting d5P c /dt = 0, based on the assumption that the CR pressure is uniform along 
magnetic field lines 5 . Then, he showed that the equilibrium state is subject to an instability, 
which is now known as the Parker instability. 

The Parker's work has since been elaborated upon. For instance, Giz & Shu (1993), 
Kim et al. (1997), and Kim & Hong (1998) investigated the modification of the Parker 
instability under nonuniform gravities. It was found that the linear growth rate increases 
under the gravities described by linear and hyperbolic-tangent functions. Kim et al. (2000) 
and Santillan et al. (2000) incorporated the multi-component nature of the ISM in a realistic 
gravity model, where the growth timescale turns out to be ~ 3 x 10 7 years and the length 
scale enlarges up to ~ 3 kpc. The effect of the irregular, random component of magnetic 
field was studied in Parker & Jopikii (2000) and Kim & Ryu (2001). With the strength of 
the random component comparable to the regular one (see, e.g., Beck et al. 1996; Zweibel & 
Heiles 1997), it was shown that the Parker instability can be completely stabilized. 

CRs form an important constituent of the ISM, with their energy density comparable to 
those of gas and magnetic field (see, e.g., Blandford & Eichler 1987). The analyses of Parker 
(1966, 1967) and Shu (1974) showed that CRs play a significant role in the development 
of the Parker instability by widening the range of unstable wavelengths and increasing the 
growth rate under the limit of k\\ — > oo (very large diffusion along magnetic field lines) and 
k±_ — (negligible diffusion across field lines). However, it is certainly true that n\\ and k± 
are finite (see §2.2 for details), but there has been no follow-up work on the effect of CRs with 
finite diffusion in the Parker instability. In an approach based on a different perspective, 
Nelson (1985) incorporated the dynamics of CRs by approximating their pressure as 

P c ,ij = P c ,±8ij - (P c ,± - Pc,\\)—^r (!) 

(for comparison, see equation (8) for the spatial diffusion tensor). Although the diffusion 
process was not included, the anisotropic nature of CR dynamics was taken into account. 
The surprising result was that the anisotropic CR pressure works towards stabilizing the 
instability, contrary to the common belief that CRs act as one of the agents to induce the 
instability itself. 



5 Later, Shu (1974) refined this treatment in a more intuitive form that the diffusion along magnetic field 
lines is very large but the diffusion across field lines is negligible. Both formulations of Parker (1966, 1967) 
and Shu (1974) turn out to be same in the linear regime. Sec §3.2. 
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In this paper, we describe an linear analysis analysis where CR dynamics is incorporated 
with the diffusion-convection equation (see, e.g., Skilling 1975), and their effect on the Parker 
instability is addressed. Previously, Kuznetsov & Ptuskin (1983) attempted a similar analy- 
sis. They argued that CRs enhance the instability, by showing that the critical value of the 
gas adiabatic index for the instability increases due CRs. Here, we first estimate the values 
of K\\ and k± which are applicable in the analysis, and then derive the dispersion relation. We 
show that the growth rate and the range of unstable wavelengths increase due to CRs. Our 
result confirms the validity of the analyses of Parker (1966) and Shu (1974) at a quantitative 
level, but disagrees with that of Nelson (1985). Recently, Hanasz & Lesch (2000) studied 
the Parker instability triggered by the CRs injected in supernova remnants. They solved 
numerically the flux tube equation for the magnetic field along with the diffusion-convection 
equation for CRs. Their work is the first example which took into account the CR diffusion 
in the Parker instability. However, it still needs to be quantified how much CR diffusion 
affects the range of unstable wavelengths and the growth rate. This paper addresses that 
specific issue. 

In §2 the stability analysis is described and the dispersion relation is driven. Also 
discussion on the CR diffusion tensor is presented. Interpretation of the dispersion relation 
is described in §3. Summary and discussion on the implications of our result are given in 
the final section. 



2. LINEAR STABILITY ANALYSIS 



2.1. Basic Equations 



The equation set for our purpose is the combination of the MHD equations and the CR 
diffusion-convection equation, 

dp 



P 



dv 
dt 



+ (v ■ V)v 



dt 



^ + v-M = o, 

f B 2 \ 1 

— = V x (v x B), 
+ v-VP g + 7g PgV • v = 0, 



dPc 
dt 



+ v-VP c + 7cP c V • v = V((«ij> VP C ) + So, 



(2) 
(3) 

(4) 
(5) 
(6) 
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where the subscripts g and c stand for gas and CRs. CRs are described by the two-fluid 
model, which is derived from the second particle momentum moment of the well-known 
CR Fokker-Planck equation (Skilling 1975). Hence, CRs are described in equation (6) by 
a pressure, plus an equation of state for the CRs represented by the adiabatic index, 7 C = 
1 + P c /E c , instead of a full momentum distribution function (see Drury & Volk 1981; Jones 
& Kang 1990, for details of the two fluid model), (/c^) is the energy weighted mean diffusion 
tensor of CRs (see §2.2 for discussion). The source term, So, in the CR pressure equation is 
introduced to set up an initial equilibrium state (see §2.3), not to describe the injection from 
thermal particles to CR particles. We ignore that process, along with shock acceleration 
of CRs (see, e.g., Drury 1983; Blandford & Eichler 1987, for details of shock acceleration). 
There are no shocks in the regime of linear stability analyses. 

As pointed out by Parker (1966, 1967), on the scale where the Parker instability is 
relevant, the dominant contribution to gas pressure would not come from the thermal motions 
of atoms or molecules, but would come from the turbulent motions of cloudlets. Then, the 
value of the "effective" adiabatic index for gas, 7 g , should be determined by considering 
the detailed mechanisms involved, such as supernova explosions, stellar winds, the Galactic 
differential rotation, cloud-cloud collisions, turbulence dissipation, and etc. Although there 
has been much progress in the studies of each mechanism, the determination of 7 g for an 
ensemble of cloudlets is less well understood. Hence, here we simply set 7 g = 1, assuming 
that the cloudlet random motions are constant. 

The adiabatic index of the CRs, 7 C = 1 + P c /E c , can be simply related to the form of 
the CR momentum distribution if the latter is a power law with an index between 4 and 5. 
In particular a momentum distribution 

f(q)(xp- q with g~14/3 (7) 

appropriate for Galactic CRs (see, e.g., Blandford & Eichler 1987), leads to P c / E c = (q — 
3)/3 ~ 5/9. Hence, 7 C = 14/9 is used in our analysis. 



2.2. Cosmic Ray Diffusion Tensor 

The frequently used form of the CR diffusion tensor is 

Kij = K,±Sij - (k± - + e ijk K A ^-, (8) 

where E>i is the magnetic field vector; k\\ and k± are the diffusion coefficients along and 
across mean field, respectively, and ha represents the curvature and gradient drifts (see, e.g., 
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Bieber & Matthaeus 1997; Giacalone & Jokipii 1999; Casse et al. 2002, for discussions on 

Kij). 

Giacalone & Jokipii (1999) and Casse et al. (2002) used Monte Carlo simulations in 
modeling turbulent magnetic fields and estimated the diffusion properties. With the energy 
ratio of the random to total magnetic fields 

SB 2 

X ~ B 2 + 5B 2 ' ( ) 



Casse et al. (2002) showed that whenever x < 1 ; Bohm diffusion (k oc p/ yp 2 + m 2 c 2 ) does 
not apply, but the quasi-linear approximation does for parallel diffusion. They found that 

H = lit ^ h - 0A X(r L krmn) 2/3 , (10) 

for the Kolmogorov turbulence. Here, r L and k min represent the Larmor radius and the 
minimum wavenumber of the Kolmogorov spectrum, respectively. The formal quasi-linear 
result would replace 0.4 by ir/Q. Note that equation (10) applies even when x — 0.99. 

Based on the GALPROP model of CR propagation in the ISM, Strong & Moskalenko 
(1998) found that k\\ ~ 6 x 10 28 cm 2 /sec at the rigidity of r L /B = 3GV using isotropic 
diffusion. Matching this with equation (10), k\\ can be estimated. Taking the momentum 
distribution in equation (7) and letting the energy weighted mean diffusion be 

/«(v/^TT-l) f(q)p 2 dp 
(«> = > f , (11) 

we get 

where A and Z are the CR atomic number and charge, respectively, and L max is the coherent 
length of the regular component of magnetic field. Note that equation (12) is derived with 
X — 0.2, but the result of Strong & Moskalenko (1998) was based on isotropic diffusion 
(x — 1)- However, the normalization seems to be uncertain, at least, by a factor of a few 
anyway, since the details of CR propagation are not really understood. 

Casse et al. (2002) found the perpendicular diffusion to behave according to 

« ± ~0.2 X 2/3 «||, (13) 

rather than k± ~ 10~ 6 /t|| which is predicted in the quasi-linear result. Also, from the 
argument based on the escape of the Galactic CRs and their lifetime, Giacalone & Jokipii 
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(1999) drew a consistent value for the perpendicular diffusion, k± ~ 0.02 — 0.04/t||. On the 
other hand, Bieber & Matthaeus (1997) argued that 

(c \ C Tl 1 
Tdecorr «± with — T dccorr ~ 2- , (14) 

where r de corr is the CR decorrelation time. So with r L /L max <C 1, it is expected that ka <C k_l 
and can be neglected to a first approximation. 

In the rest of the paper, brackets will be dropped in the mean diffusion coefficients, for 
simplicity. 



2.3. Initial Equilibrium State 

A stability analysis is started by setting up the initial equilibrium configuration. We 
employ the one originally suggested by Parker (1966, 1967). In the Cartesian coordinates 
(x, y, z), the azimuthal magnetic field is set to lie along the y-direction (0, -BoN> 0), and the 
externally given uniform gravity to accelerate in the negative z-direction (0,0,— g). Then, 
the initial state of mass density, po, gas pressure, P g o, CR pressure, P c o, and magnetic field, 
B , are described by an exponential function 



po(z) _ P g0 (z) _ P c0 (z) _ Bj(z) 
p (0) P g0 (0) P c0 (0) P 2 (0) 




where H = (1 + a + /3)a 2 /g and a is the isothermal sound speed (with 7 g = 1). The scale 
height (160 pc) and the velocity dispersion (6.4 km s _1 ) of interstellar clouds will be used 
for H and a (see, e.g., Falgarone & Lequeux 1973). a is the ratio of initial magnetic to gas 
pressures and (5 is the ratio of initial CR to gas pressures, respectively, and they are assumed 
to be constant. 

Special attention needs to be given to the initial equilibrium of equation (6). Non- 
zero k± would cause CRs to diffuse upwards. The source term, So = —k±P c o(z)/H 2 , was 
included to balance it. This ad hoc treatment, however, would introduce spurious features 
in the stability properties, and the interpretation of it should be done with caution (see §3.3 
and 3.4). Of course, with k±_ = 0, this problem does not appear. 



2.4. Linearized Perturbation Equations 



Here, we focus as in Parker (1966) on the stability in the (y,z) plane defined by the 
initial magnetic field and gravity. The analysis becomes simplified, if the linearization of 
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equations (2)-(6) is proceeded with dimensionless quantities (Shu 1974). With H and H/a 
as the normalization units of length and time, we define dimensionless coordinates and time, 

y' = y/H; z' = z/H; t' = at/H, (16) 

and introduce dimensionless perturbations of density, s, velocity, u, magnetic field, b, gas 
pressure, p g , and CR pressure, p c . Then, the perturbed state can be written as 

p = p (z)(l + s); v = au; B = B (z)(e y + b); P g = P g0 (z)(l + p g ); P c = P c0 (z)(l + p c ), 

(17) 

where e y is the unit vector along the ^/-direction. For simplicity, primes have been dropped 
in equation (17) and will be in the rest of the paper. 

Substituting equation (17) into equations (2)- (6) and keeping terms only up to the linear 
order of perturbations, the linearized perturbation equations are written as follows: 

ds du v du z . . 

^ + ^ + ^ = (18) 

OUy dp, dPc ( ) 

dt + dy +IJ dy +Mz U ' 1 9j 

+ (1 + a + P)s - (p g + [3p c + 2ab y ) + — (p g + (3p c + 2ab y ) - 2a-^ = 0, (20) 

db v 1 du z . . 

9f-2"' + a7 = ' (21) 

db z du z n 

■w~W = ' (22) 

dp g f du y : du 



m "- +1 H% + ^ 1=0 - (23) 



dp c (du y du z \ d 2 p c ( dp c d 2 p c \ db z 

Here, again primes have been dropped in normalized km and k± for simplicity. 
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2.5. Dispersion Relation 

The normal mode of perturbations has the following form 



s(y,z,t) 




s 


u y (y,z,t) 




Uy 


u z (y,z,t) 




u z 


b y (y,z,t) 




by 


b z (y,z,t) 




b z 


p g (y,z,t) 




Ps 


. p c (y,z,t) _ 




. Pc . 



exp(nt) exp(—irjy) exp I -z — i(z 



(25) 



where n is the dimensionless growth rate, r\ and ( are the dimensionless wavenumbers along 
the azimuthal (y) and vertical (z) directions, respectively. The exp(z/2) factor was included 
due to the stratified background. The same notations were used for perturbations themselves 
in the left side and their amplitudes in the right side, because no confusion arises in later 
algebra. Substituting equation (25) into equations (18)-(24) results in the following set of 
equations: 

ns - irjUy - ( - + i( ) u z = 0, 



nu,, 



ir]p g - ir]Pp c + ab z = 0, 



nu * +{i+a+ " )s - 6 + K ) (r " + * + 2ab ' ] + n ' 2ah = °- 



nb y — i(u 2 
nb z + ii]u z 



0, 

0, 



np g - irj-f g u y - 



l-7g( 



0, 



(26) 

(27) 

(28) 

(29) 
(30) 

(31) 



n + Knrj 



H 1 2 + ^ 







p c - iri^ c u y - 


_l-7,(i-<) 



u z + (k± — K\\)irjb z = 0. (32) 



The dispersion relation is derived by combining the above seven equations. Although 
straightforward, it entails tedious algebra. Here, we present a few intermediate steps. First, 
by eliminating b z in equation (32) with equation (30), we have 



n + - k ± [ - + i( 



np c - irjlcnuy 



' ' (33) 



- 9- 



which expresses p c in terms of u y and u z . Then, upon substituting s in equation (26), b y in 
equation (29), b z in equation (30), and p c in equation (33) into equations (27) and (28), we 
obtain two equations for u y and u z 



(n 2 + 7g r? 2 ) 



n 



+ K\\r) 2 - k ± [ - + i( 



+ PjcV 2 n \u y -ir) 



n 



+ k\\V -K±[- + iC 



X 



1 + a - 7 g ( - - < 



+ (3n 



l-7c 



<^ - («± - «||)/Vj = 0, 



(34) 



and 



+ ■ 



+ «||?7 - «± - + < 



1 + a + /3 



- /?7c™ ( - + K J r" M ?/ 



n 2 + 2ar7 2 + /5 ( i + < ) + (2« + 7g ) ( 7 + C 



4 



l-7c( 2"< 



n > >m z = 0. 



(35) 



Finally, by combining equations (34) and (35), we get the dispersion relation, which is a 6th 
order polynomial of n 

n 6 + C 5 n 5 + C 4 n 4 + C 3 n 3 + C 2 n 2 + C x n + C = 0, 



where 



C R = 2 



C 4 = (2a + 7g + 07c) (V + C 2 + + ^||^ 2 - kl Q + 



C 3 = (4a + 2 7g + 7c ) (V + C 2 + ^ 



kii^ 2 - «± ( - + K 



+P*± ( - + i( 



(36) 

(37) 
(38) 



(39) 



Co = 



2a( 7g + (3 lc ) (r? + C 2 + ^ - (1 + a + (5) (1 + a + (3 - 7g - /3 7c ) 
+ (2« + 7g) (V + C 2 + ^ K||77 2 -« ± Q + iC^ 



- 2+ k 



(40) 
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Ci 



x 



2a(2 7g + /3 7c ) ^ + C 2 + - \-{l + a + (3){2 + 2a + (3-2 lg -(3 lc 



rf + {3(1 + a + P)r] A (K ± - km) 



K\\V -K±[- + iC 



(41) 



2«7g ( V 2 + C 2 + I ) - (1 + a + (3){1 + a - 7g ) 



Co = 

+ /3(l + a + /3)(«L-K|i) 



K||77 - K± - +< 



(42) 



3. Results 
3.1. Parameters 

The dispersion relation, equation (36) augmented by equations (37-42), gives the linear 
growth rate, n, as a function of the azimuthal wave number, 77, and the vertical wavenumber, 
(. It involves the parameters a, (3, 7g , 7c , k\\, and k±. In addition to 7g = 1 and 7c = 14/9 
which were specified in §2.1, a — 1 and (3=1 will be used in calculation of the growth rate. 
Setting a — 1 and {3 = 1 means initially the magnetic and CR pressures are same as the 
gas pressure. For the energy weighted CR diffusion coefficients, km = 3 x 10 28 cm 2 s _1 and 
k± = 0.02k\\ will be taken as the fiducial values (§2.2). Then, in units of aH, approximately 
(k\\,k±) = (100,2). Other values of (k\\,k±) will be also considered in demonstrating the 
effect of CR diffusion. 



3.2. Parallel Diffusion (Non-zero k\\ and n± = 0) 

We first check whether our dispersion relation recovers the ones of previous works with 
simplified treatments of CR dynamics. dSP c /dt = of Parker (1966) is translated to 7c = 
and km = k± = in our formulation. The assumption introduced by Shu (1974) corresponds 
to K|| — > 00 with Kj_ = 0. Although they look different, both limits end with the same 
dispersion relation 



n 4 +(2a+ 7g )h 2 + C 2 + -K+ 



2«7g W + C + 7 " (! + « + P) (1 + « + P ~ 7g) 



which matches with those of Parker (1966) and Shu (1974). 



rf = 0, 
(43) 
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Fig. 1. — Dispersion relations of the Parker instability with non-zero k\\ and k± = 0. The 
growth rate (the largest n) is presented as a function of the wavenumber along the initial 
magnetic field direction. The vertical wavenumber along the direction of gravity was set 
to be zero. The normalization units of time and length are 2.4 x 10 7 yrs (H/a) and 160 
pc (H), respectively. Each curve is labeled by the value of k». Values of other parameters 
are specified within the frame. The growth rate of the case without CRs (/3 = 0) is also 
presented for comparison, where the same normalization was applied. 
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We now study the case of non-zero parallel diffusion, but no perpendicular diffusion. 
Note that with k± = 0, the initial equilibrium is exact, and the ad hoc inclusion of a source 
is not necessary (see §2.3). Hence, the dispersion relation is exact with its limit. Taking 
k± = 0, equations (36)-(42) reduce to 

n 5 + /t||77 2 n 4 + (2a + 7g + he) (r? + C 2 + jj n 3 + (2a + 7 g ) (V + C 2 + ^ K^n 2 



+ 
+ 



77 2 n 



2« 7 g (V + C 2 + ~ (1 + ot + (5) (1 + a + /3 - 7g ) 



k iiV 4 = 0, (44) 



after factoring out n + /t||?7 2 = which represents diffusive decay. 

Figure 1 plots the growth rate (the largest n) as a function of the azimuthal wavenumber, 
r), for zero vertical wavenumber ( = and several different values of k\\ including k\\ — > oo 
and K|| = 0. Non-zero ( reduces n. We note that the dispersion relation (44) is of fifth 
order and has five roots. As rey — > oo, four reduce to the roots of (43), while the last one 
becomes n + K\\i] 2 = representing another mode of diffusive decay. Figure 1 also plots 
the growth rate of the case without CRs (f3 = 0). The figure shows that the growth rate 
increases with increasing k». Finiteness of k» reduces the growth of the instability over that 
of km — > oo, because there is a gradient of CR pressure along magnetic field lines. The 
gradient hinders the falling motion of gas from arc regions to valleys and slows down the 
development of the instability. However, with the value of /cy expected in the ISM, ~ 100, 
the growth rate is close to that of n\\ — > oo. So we conclude that dSP c /dt — or k\\ — > oo with 
k± = used in previous analyses were good approximations and produced quantitatively 
correct results. With k\\ = 100, the maximum growth rate and the critical wavenumber 
where n — > are about twice larger than those of (5 — 0, indicating CRs can enhance the 
instability significantly. 

One thing to note is that the critical wavenumber is independent of k\\, once k\\ > 0, 
and is given as 

V 2 C = (1 + a + + a + [3 - 7g )/(2a 7g ) (45) 

The reason is the following. It takes infinitely long for such a marginally stable state to be 
developed, because of the zero growth rate. This means that, however small the diffusion is, 
the system has enough time to diffuse across any gradient of CR pressure along magnetic 
field lines. However, with no diffusion, k\\ = 0, the critical wavenumber is different and given 

as 

= (1 + a + (3){1 + a + (3 - 7g - /H)/[2a( 7g + (3 lc )] - ( 2 - -. (46) 
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which is much smaller for the parameters we employed. And as the matter of fact, the 
growth rate is much smaller too, and even smaller than that of /3 — 0. This is because 
without diffusion, CRs accumulate at valleys along with gas, and their pressure pushes gas 
out of valleys exerting a stabilizing effect. 



3.3. Perpendicular Diffusion (km = and Non-zero k±) 



We start to investigate the effect of non-zero perpendicular diffusion on the Parker 
instability, by looking at the case of k\\ — r\ — 0. This case describes the acoustic instability 
of CR mediated gas, which was studied by Drury & Falle (1986) and Kang et al. (1992), but 
with magnetic field. From the full dispersion relation (36)-(42), we get in this case 



n 



K ± 



n 



K ± 



K ± 



2 s 



n 



+ (2a + 7g + f3 lc ) 



n 



- 2 + K 



0. (47) 



In the limit of zero diffusion, the dispersion relation, n = rki-sj (2a + 7 g + /?7 c )(l/4 + ( 2 ), 
describes a pair of magnetosonic waves propagating upwards and downwards. The term of 
/37c represents the effect of additional support of CR pressure on top of the magnetosonic 
waves. The 1/4 term comes from the exp(z/2) factor in the normal mode of perturbations 
in equation (25), added to account the stratified background. As a matter of fact, all the 
1/2 and 1/4 terms in the dispersion relation (47) come from the same origin. 

Figure 2 plots the growth rate as a function of the vertical wavenumber, (, for two 
non-zero k± , s. For other parameters, the values listed in §3.1 were used. As pointed out 
in Drury & Falle (1986) and Kang et al. (1992), due to the ad hoc source term, which is 
necessary to sustain the initial equilibrium state, the analysis is justified rigorously only in 
the limit of k z ^> 1/H or ( ^> 1 in our normalized units. In that limit, the growth rate of 
non-trivial modes is 

/97c ^ P 



n 



±V2« + 7g C 



(48) 



2k ± 2 v /2a + 7 g ' 

The first term on the right hand side represents magnetosonic waves, which are shown in the 
bottom panels of Figure 2. With non-zero diffusion, the CR perturbation associated with 
with small wavelengths is wiped out due to the diffusive nature. So CR pressure (f3) does 
not contribute to the speed of the waves. The second and third terms are attributed to the 
acoustic instability, if their sum is positive, or if k± is sufficiently large, which is shown in 
the upper left panel of Figure 2. The above limiting growth rate matches exactly with that 
of Kang et al. (1992), if the magnetosonic waves are replaced by sound waves. 




Fig. 2. — Dispersion relations with k\\ = and non-zero k± for zero azimuthal wavenumber, 
77 = 0. The real part, n r , and imaginary part, rij, of the grow rate are presented as a function 
of the vertical wavenumber, (, for two different values of k± , s. The normalization units 
of time and length are 2.4 x 10 7 yrs (H/a) and 160 pc (H), respectively. Values of other 
parameters are a — 1, (3 — 1, 7 g = 1, and 7 C = 14/9. 
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The physical mechanism of the acoustic instability is the following. Perturbations in 
a gas mediated by CRs with a gradient can become unstable, because the CR pressure 
perturbation is reduced by diffusion. In the limit of large diffusion, CRs are completely 
decoupled from the gas for small scale perturbations but the gradient of CR pressure will 
remain the same. In this limit, let us suppose a constant volume force, F, is exerted on 
the gas so that the acceleration is F/p. Then, with —FSp/p 2 , compressed regions in a wave 
train will be accelerated in the opposite direction of the applied force, while decompressed 
regions will be accelerated in the direction of the force. As a result, oscillating density 
disturbances (or magnetosonic/sound waves) moving opposite to the direction of the force 
will suffer an extra restoring force and their amplitude will grow. On the other hand, waves 
propagating in the other direction will be damped. This explains why disturbances traveling 
in one direction will be amplified, while those traveling in the opposite direction will decay. 
The instability works only if the following two conditions are satisfied: 1) the perturbation 
wavelength is shorter than the scale height of CR pressure, or r] > 1, and 2) the scale height 
of CR pressure is smaller than the diffusion length associated with sound speed, or k± > 1, 
in our normalized units. 

In addition to the acoustic instability (upper left panel), Figure 2 shows the unstable 
nature in the regime of k± ~ 1 and ( < 1, especially in ( = which corresponds to the 
state of no perturbation at all (upper right panel). Without any viable mechanism, we 
attribute this to the artifact of the ad hoc source term along with the exp(z/2) factor in the 
normal mode. The detailed analysis shows that it is dominated by the perturbation of CRs, 
confirming that it is not related to the acoustic instability. 

3.4. Non-zero Parallel and Perpendicular Diffusions 

Finally, we consider the effect of non-zero perpendicular diffusion on top of non-zero 
parallel diffusion. Figure 3 plots, from the dispersion relation (36)- (42), the growth rate (the 
largest n) as a function of the azimuthal wavenumber, 77, for two different sets of the values 
of k's and zero vertical wavenumber ( = 0. Comparing Figure 1 and Figure 3, it can be 
seen that although k± is expected to reduce n, the growth rates with the same k\\ are almost 
identical except around 7] = 0. This is because k± <C k» in the ISM. The peak at ( = is 
again attributed to the artifact of the ad hoc source term along with the exp(z/2) factor in 
the normal mode, as explained in the previous section. So we conclude that with realistic 
values of CR diffusion in the ISM, the effect of non-zero perpendicular diffusion is mostly 
negligible. 
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Fig. 3. — Dispersion relations of the Parker instability with non-zero k\\ and non-zero k±. 
The growth rate (the largest n) is presented as a function of the wavenumber along the 
initial magnetic field direction. The vertical wavenumber along the direction of gravity was 
set to be zero. The normalization units of time and length are 2.4 x 10 7 yrs (H/a) and 160 
pc (H), respectively. Each curve is labeled by the values of re's. Values of other parameters 
are specified within the frame. 
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4. SUMMARY AND DISCUSSIONS 

The Parker instability in the ISM is induced by the buoyancy of magnetic field as well 
as CRs. Hence, its analysis can be completed with full treatment of CR dynamics. However, 
previous analyses incorporated CRs with simplified assumptions or ignored CRs completely, 
partially due to the lack of available CR physics. For instance, although diffusion of CRs is 
finite in the ISM, Parker (1966) and Shu (1974) assumed that the diffusion along magnetic 
field lines is large enough that there is no CR pressure gradient along them, while the 
diffusion across field lines is neglected. In this contribution, we have relaxed this assumption 
and studied the role of full CR dynamics with finite CR diffusion in the Parker instability. For 
it, first, the values of the energy weighted mean diffusion coefficients, which are applicable 
to the scales relevant to the Parker instability, have been estimated as k\\ ~ 3 x 10 28 cm 2 
s _1 and k± = 0.02aC|| . Then, a standard normal mode analysis has been performed in the 
two-dimensional plane defined by the gravity and the initial magnetic field. Linearized 
perturbation equations have been combined into the dispersion relation (36) of a polynomial 
of 6th order in the growth rate n with complex coefficients given by (37)-(42). 

It has been shown that finiteness of parallel diffusion slows down the development of the 
Parker instability. However, with k\\ — 3 x 10 28 cm 2 s _1 in the ISM, the maximum growth 
rate is smaller only by a couple of percents than that for k\\ — > oo, and the range of unstable 
wavenumbers remains the same. Inclusion of perpendicular diffusion with k,± = 0.02k\\ 
doesn't change the growth rate noticeably. That is, the original Parker's approximation 
of infinite k\\ and k± = was a good one and produced a quantitatively correct result. 
Hence, we conclude that CRs can enhance the Parker instability significantly, by increasing 
the maximum growth rate by a factor of up to two or so. We would like to note that this 
result disagrees with that of Nelson (1985), who found that CRs could stabilize, rather than 
destabilize, the Parker instability. 

As noted in the introduction, recent studies of the Parker instability, where the random 
component as well as the regular component of magnetic field were considered (Parker & 
Jopikii 2000; Kim & Ryu 2001), showed that the random component of strength comparable 
to that of regular component (5B 2 /B 2 > 0.5) can stabilize the instability completely. For 
smaller 5B 2 /B 2 , the instability is still operating, but with reduced growth rate and vanishing 
wavenumber along the radial direction of the Galaxy. Hence, it was concluded that the Parker 
instability alone has a difficult time forming Galactic structures like giant molecular clouds. 
However, our new result suggests that the Parker instability might be preserved once CRs 
are incorporated, since inclusion of CRs increases not only the growth rate and but also 
the range of unstable wavenumbers. Settling this issue would require full three-dimensional 
analysis with random component of magnetic field, which we leave for future work. 
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Finally, we comment on the applicability of an analysis that assumes a smooth distri- 
bution of CRs, even though they are thought to be products of discrete sources; namely, of 
supernova remnants. For the purposes of the present calculation the approximation should 
be quite adequate, in fact. Observations of secondary pion-produced 7-rays show that the 
galactic hadronic CR distribution is smooth on large scales (see, e.g., Bloeman et al. 1986). 
This is very reasonable in light of the long containment time of such CRs in the galaxy 
(~ 10 7 yrs; see, e.g., Connell 1998) and their associated diffusion and advection. Recent, 
sophisticated models of the CR distribution including stochastic sources along with numer- 
ous experimental constraints lead above a few hundred MeV to very smooth distributions 
outside the CR acceleration sites (Strong & Moskalenko 2001). Further, there are good ar- 
guments that small, isolated supernova remnants expanding into dense media could be far 
less common and less important sources of CRs than supernova remnants inside large, low 
density bubbles where the freshly accelerated CRs would be more broadly distributed (see, 
e.g., Higdon et al. 1998). 

The work was supported in part by KRF through grant KRF-2000-015-DS0046. We 
thank the anonymous referee for constructive comments. 
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